Magnetic Excitations of Stripes and Checkerboards in the Cuprates 
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We discuss the magnetic excitations of well-ordered stripe and checkerboard phases, including 
the high energy magnetic excitations of recent interest and possible connections to the "resonance 
peak" in cuprate superconductors. Using a suitably parametrized Heisenberg model and spin wave 
theory, we study a variety of magnetically ordered configurations, including vertical and diagonal 
site- and bond-centered stripes and simple checkerboards. We calculate the expected neutron scat- 
tering intensities as a function of energy and momentum. At zero frequency, the satellite peaks of 
even square- wave stripes are suppressed by as much as a factor of 34 below the intensity of the main 
incommensurate peaks. We further find that at low energy, spin wave cones may not always be 
resolvable experimentally. Rather, the intensity as a function of position around the cone depends 
strongly on the coupling across the stripe domain walls. At intermediate energy, we find a saddle- 
point at (-7T, 7r) for a range of couplings, and discuss its possible connection to the "resonance peak" 
observed in neutron scattering experiments on cuprate superconductors. At high energy, various 
structures are possible as a function of coupling strength and configuration, including a high energy 
square-shaped continuum originally attributed to the quantum excitations of spin ladders. On the 
other hand, we find that simple checkerboard patterns are inconsistent with experimental results 
from neutron scattering. 

PACS numbers: 74.72.-h, 74.25.Ha, 75.30.Ds, 76.50.+g 
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I. INTRODUCTION 

Strong correlations in electronic systems, and espe- 
cially competing interactions, can cause mesoscale elec- 
tronic structure to spontaneously develop. In cuprate su- 
perconductors and the related nickelate compounds, sev- 
eral locally inhomogenous electronic phases have been 
proposed, involving charge order, spin order, or both. 
Several experimental probes corroborate some level of lo- 
cal order, including STMjii 2 ^ neturon scattering^^ii, 
NMR 8 and /iSR studies^ Recent experimental advances 
have made possible the detection of high energy neutron 
scattering spectra^S*^ There has been much interest in 
the interpretation of these spectra in the cuprates, espe- 
cially since the high energy results from neutron scat- 
tering on La2- 2; Ba a ;Cu04 (LBCO) and YB^CujOe+a 
(YBCO) exhibit universal behavior^4& This paper ex- 
tends the previous work of Refsi n i 12 to look at a broader 
range of ordered structures, and also to explore the scat- 
tering patterns which are possible at high energies. It has 
been suggested, for example, that the high energy mag- 
netic excitations in LBCO and YBCO may be due to the 
quantum excitations of stripesAi&iiiSi We have reported 
elsewhere^ that high energy excitations near a quantum 
critical point (QCP) to disordered ladder o 13 i 14 i 15 i 16 can 
strongly resemble semiclassical excitations, due to the 
small critical exponent r\ = 0.037 associated with this 
QCP. 

One continuing mystery about the low energy results 
has been the lack of observed spin satellite peaks in neu- 
tron scattering, and also that spin wave cones are rarely 
observed in the cuprates. Rather, what is often seen in 
the low energy regime may be more accurately termed 



"legs of scattering". Both of these results have raised 
questions about a "stripe" interpretation of the data. 
We show below that at uj = 0, satellite peaks for even 
the most extreme case of square-wave spin stripes have 
very low intensity, and may not be resolvable without 
very high experimental resolution. In addition, although 
a spin ordered state results in spin wave cones due to 
Goldstone's theorem, the intensity is not always uniform. 
Rather, the intensity can be gathered on the inner branch 
of the spin wave cones [the side nearest (it, 7t)], or on the 
outer branch, depending upon the relative strength of the 
spin coupling across the charge stripes. For this reason, 
while spin wave cones are always present for ordered spin 
stripes, they may not yet be resolvable experimentally. 

An important point we wish to emphasize is that al- 
though stripes are a unidirectional modulation in an oth- 
erwise antiferromagnetic texture, they are a fully two- 
dimensional (2D) spin order, with a 2D magnetic Bra- 
vais lattice, which gives rise to 2D scattering signals at 
all energies. Our results never show streaks of scattering 
in the low energy structure, although streaks are possi- 
ble in the high energy structure for weak coupling, as we 
report below. 



II. MODEL AND METHOD 

In this work we concentrate on static stripes and 
checkerboards as arrays of antiphase domain walls in an 
otherwise antiferromagnetic texture. We are interested 
solely in the response of the spin degrees of freedom. The 
dynamics of the charge component which must reside on 
every domain wall is to renormalize the effective spin 
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(b) VB4 
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(a) CS3 



(b) CS4 
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(c) DS3 
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(d) DB2 



FIG. 1: (Color online) Stripe configurations studied in this 
paper. The coupling between nearest neighbor spins is J a , 
and the coupling between spins across a domain wall is Jb, 
as indicated in the figures, (a) VS4: Vertical site-centered 
stripes with spacing d = 4. (b) VB4: Vertical bond-centered 
stripes with spacing d = 4. (c) DS3: Diagonal site-centered 
stripes with spacing d — 3. (d) DB2: Diagonal bond-centered 
stripes with spacing d = 2. 



couplings in the Hamiltonian. We thus use a suitably 
parametrized Heisenberg model to describe the elemen- 
tary excitations of spin stripes, employing a modulation 
of the exchange integral to capture the effective spin cou- 
pling in a well ordered stripe or checkerboard phase. 



H -\Y1 Jr > r ' 



(1) 



where J r . r ' is the effective spin coupling. We work in 
units where h = 1. Nearest neighbor couplings are 
antiferromagnetic with J r r i = J a > within each 
antiferromagnetic patch. Couplings across a domain 
wall are different, and depend upon the configuration 
(such as spacing and direction of the domain walls). 
These are enumerated below. When comparing to, e.g., 
the cuprates (nickelates), our lattice corresponds to the 
copper (nickel) sites within the copper-oxygen (nickel- 
oxygen) planes. With these materials in mind, we restrict 
ourselves to patterns embedded in a two-dimensional 
square lattice. 



A. Stripe and Checkerboard Configurations 

In the presence of a host crystal, stripes are constrained 
by the symmetry of the crystal to run along major crys- 



7"? 

a 


n 

H 


i t 
jp 


t i 
i t 

i t 
t i 


t + 


i t 


t i 


t i 
i t 



it i fTf i t 
u t m t * 

It 1 t ! t It 
ft i""t"ft i t" 
U t i 1 t i 
i t 4 t | t 1 t 


t 1 t 

4 t 1 

t + t 

_ 

m> 4 

WJ 11 


t i t^r 

u t i ur* 

it 1 t i t l t 


t 4 t 

4 t 4 i 
t 4 t i 



(c) CB2 



(d) CB3 



FIG. 2: (Color online) Checkerboard patterns. The coupling 
between nearest neighbor spins is J a , and the coupling be- 
tween spins across a domain wall is Jb, as indicated in the 
figures, (a) CS3: Checkerboard sited-centered pattern with 
spacing d = 3. (b) CS4: Checkerboard site-centered pattern 
with spacing d = 4. (c) CB2: Checkerboard bond-centered 
pattern with spacing d — 2. (d)CB3: Checkerboard bond- 
centered pattern with spacing d — 3. 



tallographic directions. We consider two classes of stripe 
patterns, vertical and diagonal, as well as checkerboard 
patterns. These classifications concern the pattern of the 
antiphase domain walls in the antiferromagnetism. We 
use the term "vertical" to describe stripes whose domain 
walls run parallel to the Cu-0 (Ni-O) bond direction, 
and the term "diagonal" to describe stripes that run 45° 
from that direction. These stripe patterns are depicted 
in Fig. □ 

A further distinction between types of stripes and 
checkerboards depends on where the antiphase domain 
walls sit with respect to the atomic lattice. "Site- 
centered" stripes have domain walls which are centered 
on the sites of the square lattice (i.e. on the nickel or 
copper sites). These have antiferromagnetic coupling 
Jr,r> = Jb > across the domain wall. On the other 
hand, when the domain wall is situated between two 
square lattice sites (i.e. between nickel or copper sites, 
on the planar oxygens), the stripes are "bond-centered". 
In this case the coupling across the domain wall is effec- 
tively ferromagnetic (to preserve the antiphase nature of 
the domain walls), and Jb < OH 



Checkerboard patterns have been proposed to explain 
the real space structure observed in STM experiments on 
Bi 2 Sr 2 CaCu 2 8+( ^ (BSCCO) and Caa-xNa^CuOzCli 1 
(Na-CCOC). It is important to note that STM is a sur- 
face probe, and to date, the 4x4 pattern observed in 
BSCCO and Na-CCOC has not been confirmed by neu- 
tron scattering or other bulk probes in these materials. 
Likewise, the 4x4 pattern has not yet been confirmed 
via STM to be present in the lanthanum and yttrium 
compounds. Nevertheless, it has been noted that the 
length scale of the charge periodicity observed in BSCCO 
and Na-CCOC via STM is half that of the spin period- 
icity found in neutron scattering in related lanthanum- 
based and yttrium-based cuprate superconductors, and 
that therefore the two probes may be observing similar 
charge and spin modulations. Since the neutron scat- 
tering shows satellite peaks around antiferromagnetism, 
rather than a peak at the antiferromagnetic wavevector 
(-7T, 7r), any universal spin texture which is consistent with 
the proposed checkerboard pattern must also incorpo- 
rate antiphase domain walls in the corresponding spin 
texture. Representative, simple spin checkerboard con- 
figurations are shown in Fig. [21 The parameter d is the 
spacing between domain walls. The low energy peaks of 
these simple checkerboard patterns are inconsistent with 
neutron scattering datajiS^ as we will show in Sec. 
although more complicated checkerboard patterns^* may 
be consistent. 

We use the notation VSd,VBd,DSd and DBd to refer 
to vertical (V) or diagonal (D) stripes of spacing d in a 
site (S)- or bond (B)-centered configuration^ as illus- 
trated in Fig. ^ The notation CSd and CBd refers to 
checkerboards of site (S)- or bond (B)-centered domain 
walls which are d lattice sites apart. 
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FIG. 3: (Color online) VS4 at J b = 0AJ a : Constant energy 
cuts with windows 0.2Ja for twinned vertical, site-centered 
stripes of spacing d = 4 at Jj, = 0.4J a . The energy E is in 
units of J a S. 



the spin system with excitation energy u>f. We report 
single magnon excitations, and neglect possible spin- wave 
interactions since they are higher order effects. 



B. Spin-wave method 

We use linear spin- wave theory to explore the semiclas- 
sical spinwave excitations of well-ordered spin stripes, as 
modeled by Eq. ^ Ladder operators may be used to 
rewrite the Hamiltonian in term of spin wave excitations 
above the semiclassical ground state: 

H = \ E JrAs z r s z r , + \{sts;, + s;s+)]. (2) 

<r,r'> 

The spin ladder operators may be approximated as 
Holstein-Primakoff bosons, a standard procedure de- 
scribed elsewhere; 11 ' 20 ^ 1 in order to calculate the spin 
wave excitation spectrum as well as the zero-temperature 
dynamical structure factor, 

S(k,u;)=£ £ \(f\S\k)\0)f8(u;-u f ) (3) 

which is proportional to the expected neutron scattering 
intensity for single magnon excitations. Here |0) is the 
magnon vacuum state and |/) denotes the final state of 



III. RESULTS: SPECTRA OF VERTICAL 
STRIPES 

In this section, we report new results on the expected 
neutron scattering intensity at constant energy for verti- 
cal and diagonal stripes. In order to make a comparison 
with the experiments, we report constant energy plots of 
intensity in the (k x ,k y ) plane. We work in tetragonal 
units, where the k x and k y directions are oriented along 
the Cu-0 (Ni-O) bond direction and the antiferromag- 
netic wavevector is at Qaf = ( 7r i 7r )- I n each plot, we 
integrate over an energy window of ±0.2 J a S. For verti- 
cal stripes, we show each plot for a twinned pattern of 
stripes, adding the intensity from domains rotated 90° 
with respect to each other. 

We first consider vertical stripes. Figs. 03 and 0] show 
results for vertical site-centered stripes at coupling ratio 
•hi J a = 0.4 and Jb/J a — 1-0, respectively. Figs. ISJandEI 
show results for vertical bond-centered stripes at coupling- 
ratio \Jb/J a \ — 0.4 and \Jb/J a \ — 1-5, respectively. 

The main incommensurate (IC) peaks indicating stripe 
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FIG. 4: (Color online) VS4 at Jb — J a : Constant energy cuts 
with windows 0.2Ja for twinned vertical, site-centered stripes 
of spacing d — 4 at Jb — J a - The energy E is in units of J a S. 
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FIG. 5: (Color online) VB4 at J b = —0AJ a : Constant energy 
cuts with windows 0.2Ja for twinned vertical, bond-centered 
stripes of spacing d = 4 at Jb = — 0.4J a . The energy E is in 
units of J a S. 



order have high intensity at low energy in Figs. 13161 How- 
ever, the low energy satellite peaks are so weak as to be 
virtually unresolvable in the site-centered case (Figs. |3] 
andQJ, while they are visible but still weak in the bond- 
centered case (Figs . [5] and [SJl - At u = 0, the ratio between 
the main incommensurate peaks and the satellite peaks 
is 34 for VS4, and it is 5.9 for VB4. This is the maximum 
intensity ratio, for, e.g., a square wave pattern that fol- 
lows Figs. H (a) and (b), in which every occupied site has 
the same spin moment. With such a dramatic suppres- 
sion of the satellite peaks for VS4, and given the fact that 
the spin wave cone emanating from the satellite peaks in 
the site-centered case rapidly diminishes in intensity with 
increasing energy^ii it is not surprising that the satellite 
cones are too faint to be seen in the lowest energies of 
Figs. 01 and 0] Rather than evidence against stripes, the 
absence of observed satellite peaks may simply be due 
to insufficient experimental resolution, especially in the 
case of site-centered stripes. Any envelope softer than 
a square wave pattern will diminish the satellite peaks 
even further^ 

At slightly higher energy (the middle left panel in 
Figs. l.'iHH . notice that the intensity on the spin wave 
cones can be different on the "inner branch," the side 
closest to (tt, 7t), and the "outer branch," the side farthest 
from (it, it). In fact, in Fig. (VB4 at J b /J a = -0.4), 
the intensity is strongest on the inner branch. On the 



other hand, at higher coupling ratios (Figs.^andEJ), the 
strongest intensity is on the outer branch. The inten- 
sity ratio between inner and outer branches is a function 
of the coupling ratio \Jb/J a \, as shown in Fig. |7| To 
obtain the figure, we have taken a cut through /c-space 
perpendicular to the stripe direction, along (k x , tt) , and 
restricted the plot to a low energy, E = 0AJ a S . We 
divide the peak intensity at the inner branch of the spin 
wave cone by the peak intensity of the outer branch of 
the spin wave cone to derive the inner to outer branch 
intensity ratio in Fig. Based on a fit of the results in 
Fig. we find the following functional form for the in- 
tensity ratio of inner to outer branches for vertical stripes 
of spacing d — 4: 



S(k in ,E ) 
S(kout,E ) 



= a + 



\Jb/Ja\ 



(4) 



where k; n denotes the wavevector of the inner branch of 
the spin wave cone on the (tt, tt) side, and k out denotes 
the wavevector of the outer branch away from (tt, tt) at 
a particular energy E a . For VS4 (site-centered stripes) 
at E = 0AJ a S wc find that a = 0.85, b = 0.076, 
c = 0, and /i = 0.89. For VB4 (bond-centered stripes) at 
E = 0AJ a S we find that a = 0.84, b = 0.18, c = -0.033, 
and n = 0.84. For both site- and bond-centered stripes, 
for small enough \Jb/J a \, the intensity ratio is so dra- 
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FIG. 6: (Color online) VB4 at J;, = — 1.5J a : Constant energy 
cuts with windows 0.2Ja for twinned vertical, bond-centered 
stripes of spacing d = 4 at Jb = — 1.5J a . The energy E is in 
units of J a S. 



matic that without sufficiently high resolution, the en- 
tire spin wave cone will not be resolvable, and instead 
only "legs of scattering" will be visible. Several neutron 
scattering experiments on the cuprates report this kind 
of behavior and our calculations here indicate that 
the behavior is simply due to weak effective spin coupling 
across the charge stripes. 

For a range of coupling ratios, the acoustic band has a 
saddlepoint at (ir, 7r) which has many similarities to the 
"resonance peak" observed in the cuprates. These sad- 
dlepoints can be seen in the bottom right panel of Figs. [3} 
Here, the intensity is gathered into one main peak at 
the antiferromagnetic wavevector Qaf, where the max- 
imum intensity is higher than that at nearby energies. 
The saddlepoint structure gives rise to an hourglass shape 
emanating from the resonance peak for twinned stripes 
in energy vs wavevector plots ^^^^2^ The resonance 
peak associated with a saddlepoint is more pronounced 
at low coupling ratio. For stronger coupling ratio, the 
weight in the spin wave cones is shifted to the outer 
branch, and the extra intensity due to the saddlepoint 
at (7r, 7r) is reduced. 

For the case of vertical, site-centered stripes of spac- 
ing p — 4 (VS4), for J b = 0.05J a , the resonance en- 
ergy is E IeB — 0.63J a 5, and for J b = 0.2J a , the sad- 
dlepoint is at _E rcs = 1.2 J a S. Assuming the value of 
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FIG. 7: Intensity ratio between the "inner" and "outer" 
branches of the spin wave cones as a function of coupling 
ratio \Jb/J a \ for vertical stripes. The ratio is calculated at 
constant low energy E = 0AJ a S. The solid points are our 
numerical calculations of the intensity ratio. The solid line is 
a fit to the results, as explained in the text, (a) VS4: Verti- 
cal, site-centered stripes of spacing d = 4. (b) VB4: Vertical, 
bond-centered stripes of spacing d — 4. 



J a (~ 140meV— ) is relatively unchanged upon doping, 
this gives a range of E les = 44 — 84meV. This encom- 
passes the range E Tes — 50 — 60meV given in Ref£ for 
the resonance peak in LBCO. For the case of vertical, 
bond-centered stripes of the same spacing (VB4), with 
Jb — 0.05J a , the saddlepoint is at E Ies = 0.6J a S, and for 
Jb = 0.2J a , £" rcs = l.U a 5. This corresponds to a range 
of E les = 42 — 77meV, again encompassing the experi- 
mental range E les = 50 — 60meV in Refi£. This would 
indicate that to match the resonance energy in LBCO 
requires that 0.05J a < J b < 0.2J a il 

Although doping is not explicitly in our model (rather 
we treat the stripe spacing in a phcnomcnological man- 
ner), low energy neutron scattering indicates that the 
stripes move closer together as doping is increasedi 28 i 29 
The energy of the saddlepoint resonance increases mono- 
tonically as the stripes are moved closer together in our 
model, holding J a and Jb fixed. If J a and Jb are rel- 
atively independent of doping (and we believe this is 



physically reasonable), then we would expect doping to 
increase the energy scale of, e.g., the resonance, as is ob- 
served in underdoped cuprates. However, predicting how 
J a and Jb depend on doping is beyond the scope of our 
present model. Another important piece of physics in 
the cuprates is that the intensity of the "resonance peak" 
increases below T c . However, since we do not explicitly 
include superconductivity, our model does not address 
the observed increase in intensity of the resonance peak 
as superconductivity onsets. 

It is also possible to tune the coupling ratio \Jb/J a \ so 
that the spin wave bands cross instead of exhibiting a 
saddlepoint. This happens at Jb/J a = 1 for site-centered 
stripes, and at \Jb/J a \ — 0.56 for bond-centered stripes. 
It is unlikely that such a crossing would be observed 
experimentally, 3 - because of the fine-tuning it would re- 
quire. The high energy structure of the spin wave re- 
sponse therefore depends on whether the coupling ratio 
is above or below these special points. For small coupling 
ratios, the generic behavior is that there is a high energy 
square-shaped continuum above the saddlepoint, with 
vertices rotated 45° from the low energy incommensurate 
peaks^ reminiscent of the universal high energy behav- 
ior recently reported in LBCO 5 and YBCO£ For very 
weak coupling, the highly anisotropic spin wave structure 
yields flat square edges in the high energy structure. As 
the coupling ratio \Jb/J a \ is increased, the square edges 
distort. For higher coupling ratios, at and beyond the 
"crossing point", the high energy cross sections display 
a variety of patterns, including circulai^i and square- 
shaped continua, as well as more complicated patterns. 
For large coupling ratios, the high energy square-shaped 
continuum has vertices that are along the same direction 
as the low energy peaks, as in Fig. 



IV. RESULTS: SPECTRA OF DIAGONAL 
STRIPES 

We now consider diagonal stripes. Diagonal stripes 
have been observed in the nickelatesj^SiS and at low dop- 
ing in the cuprates, e.g., for x < 0.05 in L^-xSr^CuC^ 
(LSCQ)i 34 i 35 When the crystals are detwinned, the di- 
agonal stripe incommensurate peaks are also untwinned 
in the cupratesi^^ To make a comparison with the ex- 
periments on diagonal stripes, we report results for un- 
twinned stripes as well as twinned stripes. 

The expected neutron scattering intensities for diago- 
nal stripes depend starkly on whether the stripes have an 
even or odd spacing. Diagonal bond-centered stripes of 
odd spacing generally display net ferromagnetism. This 
is because diagonal bond-centered domain walls have a 
net magnetic moment, and at odd domain wall spac- 
ing, each domain wall will have the same moment. This 
dramatically changes the nature of the Goldstone modes 
from linear to quadratic in \k — k \. For site- or bond- 
centered diagonal stripes of even spacing, the number 
of magnetic reciprocal lattice vectors doubles, and then 
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FIG. 8: (Color online) DB2 at J b = -0.1J a : Constant en- 
ergy cuts with windows of 0.2Ja for diagonal, bond-centered 
stripes of spacing d — 2 with Jb = — 0.1J a - The left column 
is untwinned and the right one is twinned. The energy E is 
in units of J a S. 

even Qaf = (tt, tt) is a reciprocal lattice vector. Since 
stripes are arrays of antiphase domain walls in the an- 
tiferromagnetic texture, there can be no net antiferro- 
magnetism, and hence no zero frequency weight at Qaf- 
Nevertheless, a spin wave cone emanates from the (tt, tt) 
point when the diagonal spacing is even. The cone has no 
weight at zero frequency, and gains a faint appearance as 
energy is increased^ Because of this unique band struc- 
ture, even-spaced diagonal stripes cannot display a sad- 
dlepoint in the acoustic band at (tt, tt). There is, however, 
a saddlepoint in the next optical band. 

We show results for diagonal stripes in Figs. 1811 II As 
with vertical stripes, the intensity profile of the spin wave 
cones at low energy is a function of | Jb / J a \ ■ For weak cou- 
pling Jb, the weight is strongly gathered near the (tt, tt) 
point. In Fig. 1111 we plot the intensity ratio between the 
inner and outer branches of the spin wave cones at small 
energy, for both DB2 (diagonal bond-centered stripes of 
spacing d — 2, Fig. II lb) and DS3 (diagonal site-centered 
stripes of spacing d — 3, Fig.lTTbh To obtain the figures, 
we have taken a cut through fc-space perpendicular to 
the stripe direction, along (k x , —k x ), and we restrict the 
plots to a low energy, E = 0AJ a S. We plot the peak in- 
tensity at the inner branch of the spin wave cone (toward 
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FIG. 9: (Color online) DS3 at J;, = 0.1J a : Constant energy 
cuts with windows of 0.2 Ja for diagonal, site-centered stripes 
of spacing d = 3 at Jb = 0.1J a . The left column is untwinned 
and the right one is twinned. The energy E is in units of J a S. 



(•7T, tt)) divided by the peak intensity of the outer branch 
of the spin wave cone (away from (tt, it)). By fitting the 
results in Fig. El f° r diagonal stripes we find that it is 
well described by the same functional form that we used 
for vertical stripes, although the fitted constants are dif- 
ferent. Based on a fit to Eqn. 0] we find that for DB2 
at E = 0AJ a S (Fig. Ek), a = 1.8, b = 0.022, c = 0, 
and \i = 2.6. There is a slight deviation from the fit near 
J b = -0.25 J a . For DS3 at E = J a S (Fig. Eh), we find 
that a = 0.74, b = 0.15, c = -0.11, and /j, = 0.61. As 
with vertical stripes, for diagonal stripes of both the site- 
and bond-centered type, at small enough coupling ratio 
\Jb/ Ja\, the intensity on the outer branch is so weak that 
only part of the spin wave cone will be visible without 
sufficient experimental resolution, indicating that diago- 
nal stripes can also display "legs of scattering" when the 
spin coupling across the charge stripe is weak. (Simi- 



FIG. 10: (Color online) DS3 at Jj, = 0.5J a : Constant energy 
cuts with windows of 0.2Ja for diagonal, site-centered stripes 
of spacing d = 3 at Jb = 0.5 J a - Each panel is twinned. The 
energy E is in units of J a S. 

larly, large coupling across the charge stripe can display 
outwardly dispersing "legs of scattering".) 

In Fig. [S] we show constant energy cuts for DB2, a 
bond-centered diagonal stripe configuration with spacing 
d = 2 between domain walls, at weak coupling across the 
charge stripes, Jb = — 0.1J a . We report untwinned inten- 
sity plots in the left column, and twinned intensity plots 
in the right column. At low energies, the twinned inten- 
sity plots show four spin wave cones dispersing from the 
incommensurate peaks, with weight concentrated near 
the (tt, tt) point . 

At higher energies in Fig. |SJ because the stripe spac- 
ing d is even, the acoustic band cannot support a sad- 
dlepoint at (tt,tt), as discussed above. Rather, it is 
the optical band that has a saddlepoint, appearing at 
E = 1.4J a S. Like the resonance peaks in vertical stripes, 
the saddlepoint^ has higher intensity, with low energy 
and high energy branches emanating from it. Due to 
the weak coupling across charge stripes |J&| <C J a , the 
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FIG. 11: Intensity ratio between the "inner" and "outer" 
branches of the spin wave cones as a function of coupling 
ratio \Jb/Ja\ for diagonal stripes. The solid points are our 
numerical calculations of the intensity ratio. The solid line 
is a fit to the results, as explained in the text, (a) DB2: 
Intensity ratio for diagonal bond-centered stripes of spacing 
d = 2 at E — OAJaS. (b) DS3: Intensity ratio for diagonal 
site-centered stripes of spacing d = 3 at E — J a S. 



45° from the low energy peaks, a familiar pattern at high 
energy. 

In Fig. we show DS3 at higher coupling ratio, 
Jb = 0.5 J a . This pattern and coupling has been shown 
to capture many essential features^ - of the data for the 
nickelate compound L^-zSr^NiO^ at x = 1/3, and to 
some extent x = 0.275 as we\lm2> To more fully compare 
with this experiment, we show results only for twinned 
stripes. Our results may be compared with Figs. 2, 5, 
and 10 of Rcf. 32 . (However note that the axes of the 
plots in Refi2£ are rotated 45° from our plots.) At low 
energy (up to E = J a S), the spin wave cones have in- 
tensity peaked on the outer branch. However, at higher 
energy, E — 2J a S, the intensity has shifted, and is now 
peaked on the inner branch. The spin wave cones remain 
remarkably circular through this energy range, despite 
the anisotropic coupling ratio. As the spin wave cones 
touch at E = 2.3J a S, there are four regions of high scat- 
tering radiating out from (n, tt). As the cones merge, 
they form a central peak at (tt, 7r), surrounded by four 
smaller regions of high scattering, rotated 45° from the 
low energy peaks. These four smaller peaks are mainly 
due to the addition of two twinned stripe patterns. At 
higher energy, these four peaks move outward (occupy- 
ing the corners of the graph in Fig. 5(f) of Refi^), while 
the peak at (jr, it) diminishes. Finally at higher energy 
E = 3.lJ a S, the peak at (tt, tt) is no longer visible. This 
indicates that rather than being a saddlepoint, as hap- 
pened for Jt <C J a , by the time the coupling ratio has 
reached Jb = 0.5J a , there is no longer a saddlepoint at 
(it, 7t). Rather, there is simply a band edge at E = 3J a S. 
This is also consistent with the fact that the (tt, tt) peak 
in this case does not display significantly higher intensity 
compared to that of nearby energy cuts. 



V. RESULTS: SPECTRA OF 
CHECKERBOARDS 



high energy branches emanating from the saddlepoint 
are rather flat, giving rise to a high energy square shaped 
continuum in the twinned plots, rotated 45° from the low 
energy peaks. 

In Fig. El we show results for DS3, diagonal site- 
centered stripes of spacing d = 3, with Jb = 0.1J a . Be- 
cause the coupling Jb is weak, the spin wave cones have 
weight gathered near (7r,7r). The low energy dispersion 
for DS3 at this coupling is much steeper than that for 
DB2 at the same coupling strength, because unlike DB2, 
DS3 has no magnetic reciprocal lattice vectors at (7r,7r), 
and so the spin wave cones of DS3 continue up in en- 
ergy to a saddlepoint at (tt, 7t). The saddlepoint occurs 
at E — 1.3 J a S, and it has extra intensity when twinned. 
As with other saddlepoints at low Jb, there are high en- 
ergy and low energy branches emanating from it, result- 
ing in an hourglass shape in E vs. k for twinned stripes. 
The weak coupling Jb <C J a gives rise to a square shaped 
continuum above the saddlepoint, with vertices rotated 



We show in Fig. El the expected scattering response 
from a typical simple checkerboard pattern. We show 
results for a site-centered checkerboard pattern of spac- 
ing d = 3, the real space pattern of which is shown in 
Fig- Et&)- At low energies, simple checkerboards which 
use the charged lines as antiphase domain walls in the 
spin texture yield zero frequency peaks (and therefore low 
energy spin wave cones) which are in the wrong direction. 
That is to say, for vertically placed domain walls (as re- 
quired by both the charge incommensurate peaks in neu- 
tron scattering on LNSCO at dopings x > 0.05^LiS and 
the STM Fourier transform peaks in BSCCOi2£) the 
spin incommensurate peaks are rotated 45° from those 
observed experimentally. 16 Simple checkerboard pattens 
like these always give IC spin peaks rotated 45° from the 
direction of the charge IC peaksji§ii& contrary to what 
has been seen experimentally. However, more compli- 
cated checkerboard patterns^ may be capable of fitting 
the low energy data, and we plan to explore the finite 
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FIG. 12: (Color online) CS3 at J 6 = 0.1J o : Constant energy 
cuts for site-centered checkerboard pattern with stripe spacing 
d — 3 at Jt — 0.1J a . The energy E is in units of J a S. 



frequency response of these in a future publication^ 

Simple bond-centered checkerboards of odd spacing 
(d G odd) suffer yet another ill: they support a net mag- 
netization as shown in Fig.^Jb). As a result, they display 
a ferromagnetic spin peak at Qp = (0, 0), and spin waves 
which are quadratic rather than linear in \k — k a \, where 
k Q is a magnetic reciprocal lattice vector. Net ferromag- 
netism is not attainable for simple vertical stripes, but 
it is possible with bond-centered diagonal stripes of odd 
spacing, as discussed in Sec. IIVI 

The spin wave cones at low energy in Fig. El have 
weight gathered on the inner branches on the side nearest 
(-7T, 7r). The intensity in these cones merges into a square- 
like pattern as energy is increased, before the band ends 
with a peak at E = lAJ a S and (it, it). The high energy 
part of the acoustic band also has incommensurate peaks 
which are in the correct direction for the low energy IC 
spin peaks, but are overwhelmed in intensity by the cen- 
tral peak at Qaf- The high energy peak at Qaf which 
marks the top of the acoustic band is unlike the reso- 
nance peak observed in the experiment, since there is no 
scattering signal emanating from it at higher energy^ 
Rather, above E — lAJ a S in Fig. ^1 there is a spin 
wave gap to a rather flat optical band. 

The form of the energy in this checkerboard configu- 
ration (CS3) may be calculated analytically. There are 
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FIG. 13: CS3 at Jb = 0.1 J a : Dispersion and intensities for a 
site-centered checkerboard pattern with stripe spacing d = 3 
at Jb = O.lJa. The upper pannel is along (k x ,n) direc- 
tion ; the lower pannel corresponds to the diagonal direction 
{k x , k x ). The energy E is in units of J a S. 



four spins in the unit cell, as shown in Fig. Ufa), and 
only two bands. Along the direction diagonal to the do- 
main walls (k x , k x ), the dispersions for the accoustic and 
optical bands are 



k-Wc . /Ti . 3k x 



(5) 
(0) 



with A = Jb/ J a - Notice that the optical band is flat in 
this direction. Along the direction parallel to the domain 
walls (fc x ,0), the dispersions are 



(t^)V2 = (l + 3A + A 2 -Acos3fc x ) (7) 



± (l + A)Vl + A 2 + 2Acos3/c 3; , 

where the — sign refers to the acoustic band, and the + 
sign refers to the optical band. The spin wave velocities 
in the diagonal and parallel directions are 



; ll 



3VA 

3VA 
2^2 



VAF, 



VAF 



(8) 
(9) 



These band structures (with numerically calculated in- 
tensities) are shown in Fig. 1131 in the (k x , it) and (k x , k x ) 
directions. There are always two bands present, although 
one is often quite weak compared to the other. The cou- 
pling at which the two bands touch each other is Jb = J a - 
Fig. EH a ) shows the gap that is present between the 
acoustic and optical bands, and that there is no high en- 
ergy structure emanating from the peak at (jr, 7r) which 
terminates the acoustic band. 
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VI. CONCLUSIONS 

In conclusion, we have studied the expected inelastic 
neutron scattering intensities from spin waves for a vari- 
ety of spin-ordered phases, including vertical and diago- 
nal stripes and a simple checkerboard pattern. We find 
that the inelastic response is very sensitive to the cou- 
pling across domain walls throughout the energy ranges 
studied. In addition, we wish to emphasize that the elas- 
tic response can also hold surprises. For example, we find 
for vertical stripes of spacing d = 4 that when the do- 
main walls are site-centered (VS4) , the ratio between the 
main IC peaks and the next satellite peak is at least 34. 
This number is independent of coupling strength, and it 
holds for the most sharply spin-ordered case of a square- 
wave profile to the antiphase domain walls. Stripes in 
a real material are expected to have a softer profile due 
to quantum effects neglected here, and this will dimin- 
ish the satellite peaks even more. For the bond-centered 
d — 4 case (VB4), the satellite peaks are easier to detect, 
with a ratio of at least 5.9 between the main IC peaks 
and the next satellite peaks, even for the most extreme 
case of a square-wave profile. However, this is still below 
the noise for current experimental resolution. 

For vertical stripes at finite but low energies, the in- 
tensity around the spin wave cones is rarely uniform. 
Rather, it is shifted either toward or away from the (tt, tt) 
point by an amount which depends on the coupling ratio 
| Jb/ J a \- For weak coupling across the charge stripes, e.g. 
\Jb\ < J a, the weight is gathered on the inner branch to- 
ward the (tt, tt) peak, indicating that weak stripe coupling 
may be responsible for the observed "legs" of scattering 
in stripe-ordered cuprates at low frequency. For stronger 
coupling across the charge stripes, e.g. |J&| > J a , the in- 
tensity shifts to the outer branch. In either case, without 
sufficient experimental resolution, only part of the spin 
wave cones will be visible at low energy. 

At intermediate energies and for weak coupling, the 
acoustic band displays a saddlepoint whose intensity pro- 
file mimics that of the "resonance peak" observed in 
cuprate superconductors. For twinned stripes, a saddle- 
point displays the characteristic hourglass shape in E 
vs. k plots seen in some experiments, with both low 
and high energy legs emanating from a resonance-like 
peak, 4 i 23 i 24 i 25 i 26 

At high energies and weak coupling Jb, vertical stripes 
show a square-shaped continuum which is rotated 45° 
from the low energy IC peaks. As we have shown 
previously^ very weak coupling captures the spin wave 
excitations of stripe-ordered LBCO at all measured en- 
ergies, and at higher energies it bears a striking resem- 
blance to the high energy excitations of YBCO. The gap 
to spin S — 1 excitations in YBCO is likely due to the 
stripes being too weakly coupled, and therefore quan- 
tum disordered. However, at high energy, the quantum 
critical excitationsSiiii 4 * 1 ^ strongly resemble the spin 
waves studied here, which may be due to the proxim- 
ity of a QCP with small critical exponent r\ = 0.037. 12 



Near the QCP, one way to distinguish semiclassical spin 
waves from quantum critical excitations is through the 
lineshapes. Whereas semiclassical spin waves produce a 
Lorentzian lineshapc, quantum criticality yields a power 
law cusp. For intermediate and larger couplings Jb, the 
high energy response can have a variety of shapes. We 
have shown here that circular continua are possible, as 
well as square-shaped continua which either have the 
same orientation as the low energy peaks, or are rotated 
45° from that direction. 

For diagonal stripes, as with vertical stripes, the low 
energy spin wave cones have intensity profiles which de- 
pend on the strength of the coupling between the stripes. 
For weak coupling, the intensity is peaked on the inner 
branches, while for strong coupling, it is peaked on the 
outer branches. In either extreme case, the entire spin- 
wave cone may not be visible in an experiment without 
sufficient resolution. At intermediate energy, weakly cou- 
pled diagonal stripes also have a saddlepoint in the acous- 
tic band, with properties akin to the resonance peak. At 
high energy, weakly coupled diagonal stripes display a 
high energy square-shaped continuum, rotated 45° from 
the low energy peaks. At higher coupling ratio \Jb/J a \, 
the saddlepoint in the acoustic band at (tt, tt) is lost, and 
a variety of high energy scattering patterns are possible. 
Diagonal bond-centered stripes of odd spacing display net 
ferromagnetism, changing the nature of Goldstone modes 
from linear to quadratic in \k — k a \. Both site- and bond- 
centered stripes of even spacing have magnetic reciprocal 
lattice vectors which include Qaf = (t, tt)- Because of 
the antiphase nature of the domain walls, zero frequency 
is forbidden at this peak, but the (very low intensity) spin 
wave cone emanating from (tt, tt) precludes a saddlepoint 
at Qaf in the acoustic band. 

For intermediate coupling Jb = 0.5J a , the expected 
scattering intensity from diagonal, site-centered stripes of 
spacing d = 3 strongly resembles that in L^-xSr^NiC^ 
at a; = 1/3 and x — 0.275, as pointed out in Rcf. 32 . 
We find for this configuration that while the intensity 
is peaked on the outer branch of the spin wave cone at 
low energy, it moves toward the inner branch as energy 
is increased. Furthermore, the spin wave cones are re- 
markably circular, despite the anisotropic coupling ratio. 
With increasing energy, the spin wave cones merge, even- 
tually gathering weight at the central peak (tt,tt). How- 
ever, rather than a saddlepoint, this coupling has a band 
edge at (tt, tt). 

Simple checkerboards, on the other hand, have low en- 
ergy spin IC peaks rotated 45° from the observed spin 
peaks in neutron scattering^ In addition, for weak cou- 
pling across the domain walls |«7&| <C J a , rather than 
showing a resonance-like saddlepoint, the acoustic band 
has a band edge at (tt,tt)^ and therefore no branches 
emanating from it at high energy. Although the simple 
checkerboard studied here is incompatible with the ex- 
perimental data, this does not rule out the possibility of 
more complicated checkerboard patterns^ which may be 
able to capture the correct orientation of the low energy 
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spin and charge IC peaks. 
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